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It is proven that, for any soft potential <p(r) characterized by a finite Fourier transform <p(k), 
the virial and energy thermodynamic routes are equivalent if the Fourier transform of the total 
correlation function divided by the density p, h(k)/p, is an arbitrary function of pj3ip(k), where (3 is 
the inverse temperature. This class of approximations includes the mean spherical approximation 
as a particular case. 



Standard statistical-mechanical formulas relate the 
thermodynamic properties of a fluid to the two-body in- 
teraction potential (p(r) and the associated radial distri- 
bution function g(r; p, f3)r- where p is the number density 
and (3 = l/fc^T is the inverse temperature. When an ap- 
proximation is used to get g(r\ p, (3), there is in general no 
guarantee that those different formulas or routes are ther- 
modynamically self-consistent. In fact, some liquid state 
theories contain one or more adjustable state-dependent 
parameters which are tuned to achieve thermodynamic 
consistency between two or more routes. This is 
the case, for instance, of the modified hypernetted- 
chain closure^ the Rogers- Young closure^ the Zerah- 
Hansen closure,— the self-consistent Ornstein-Zernike 
approximation, 5 the hierarchical reference theory^ Lee's 
theory based on the zero-separation theorems,— the gen- 
eralized mean spherical approximation^ or the rational- 
function approximation.— It is then remarkable when an 
approximate theory for g(r; p, (3) satisfies a condition 
of thermodynamic consistency without being forced to 
do so. In this context, it is perhaps not sufficiently 
well known that the hypernetted-chain (HNC) integral 
equation provides thermodynamically consistent results 
through the virial and energy routes, regardless of the po- 
tential ip(r)J£ More recently, Mladek et al^ have shown 
that the mean spherical approximation (MSA) is exactly 
solvable for the Gaussian core model and have found 
that the virial and energy routes to thermodynamics are 
equivalent in that case as well. 

The aim of this Note is to place Mladek et al.'s 
findingii in a broader context by proving the thermody- 
namic consistency between the virial and energy routes 
for (a) any "soft" potential and (b) within a class of ap- 
proximations that includes the MSA as a particular case. 

The virial and energy routes to thermodynamics read^ 



g(r;p,fi)r ■ Vtp(r), (1) 



u(p,0) 



d 



p / dr g(r; p, 0)<p(r) 
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respectively, where p is the pressure, Z is the compress- 
ibility factor, d is the dimensionality of the system, and 
u is the internal energy per particle. The condition of 



thermodynamic consistency between both routes is 



(3) 



Let us now consider an interaction potential verifying 
the boundary conditions 

lim r d ip(r) = 0, lim r d <p(r) = 0. (4) 

r — >0 r — >oo 

While the second condition means that the potential is 
sufficiently short ranged, the first condition defines the 
kind of soft potentials to be considered here. It includes 
bounded potentials (such as the Gaussian core models or 
the penetrable sphere model^) , logarithmically diverging 
potentials jA* or even potentials diverging algebraically as 
tp(r) ~ r~ n with n < d. On the other hand, conven- 
tional molecular models (such as hard spheres, square- 
well fluids, and Lennard- Jones fluids) are excluded from 
the class of potentials ([3]). Equation Q implies that 
£>(0) = finite, where the Fourier transform of the poten- 
tial is 

$(k) = F[<p(r)\ = J dre-' k Xr). (5) 

Introducing the total correlation function h(r; p, (3) 
i i [ r : p,0) — 1, Eqs. ([T]) and ||2J| can be rewritten as 



Z(j>,0) 



■-p0<p(Q)-—p0I v {p,0), 



(6) 



u{p^) = ^+ l - P m + \ph(P,f3), (7) 

where 

I v {p,(3) = J drh{r;p,(3)r-V^(r) 

= -j^y J dkh(k;p,/3)-^-[k$(k)], (8) 

I e (p,(3) = J drh(r;p,p)<p(r) 

= tAs / dk*(*; (9) 



In the last equalities of Eqs. JSJ) and ([9]), h(k\p,(3) — 
!F[h(r\ p, (3)] and standard steps have been followed. The 
consistency condition ((3} becomes 

A(p,(3) ee A [p/ e (p,/?)] + I A = 0. (10) 



Equation (| 10[) is of course satisfied if the exact 
h(k;p,/3) is used to evaluate I v (p,(3) and I e (p,f3). On 
the other hand, as proven below, the same happens with 
any approximate function h(k\ p, f3) which, when divided 
by p, depends on p, (3, and k through the scaled variable 
z = pf3tp(k) only, i.e., 



h(k;p,p)=p- 1 F(z), z = p(3v(k), 



(11) 



where the function F(z) does not need to be specified. 
The scaling form (fTTj) implies 



d_ 

df3 



>M/,:,, .') = h(k;p,/3) + — [ph(k;p,0) . (12) 



d_ 

dp 



ph(k; p, /?)] k • _£(fc) = p(k)k • — fc(fc; p, /?). (13) 



Applying Eq. $3$ in the definition of A(p,/3), Eq. PU|) . 
one gets 



1 



d(27r) 



dk f/d>; P ,/?)A-[k^)] 



ph(k;p, (3) 



(14) 



Finally, use of Eq. (fT5|) yields A(p, /3) = 0, what proves 
the thermodynamic consistency between Eqs. ([T]) and @ 
within the class of approximations (jlip. Furthermore, Z 
and u adopt the forms 

Z(p,/3) = l+^(0)a-/3C.(a), (15) 



pu(p,p) = l + ±<p(O)a-0Ce(a), (16) 



where a ee ^(a) ee pl v (p, /3)/2d, and Ce(a) = 
— pl e (p, (3)/2. Equation (flT)]) implies the relation 
a^(a) = d [aCi,(o;)] /da, where Ce( a ) = dCe(a)/da. 

It must be noted that the compressibility route is in- 
consistent with the virial and energy routes, i.e., [1 + 
ph(0; p, /3)] _1 ^ d(pZ)/dp, for approximations of the 
form (fTTj) . While the (reduced) isothermal compressibil- 
ity 1 + ph(0; p,(3) = 1 + F ((p(Q)a) depends on p and 
(3 through the product a = pf3 only, Eq. (fl~5f yields 
d(pZ)/dp = l+cp(0)a— (3aC,' e {a), the last term depending 
on both a and (3. 

The MSA for soft potentials-^ consists of assuming 
the random-phase approximation c(r; p, (3) = —f3ip{r) 
for the direct correlation function at any distance. The 
Ornstcin-Zcrnike relation^ then gives Eq. (|11[) with 
F(z) = —z/(l + z), so that this version of the MSA is 
thermodynamically consistent through the virial and en- 
ergy routes. However, they are inconsistent with the the 
compressibility route: [1 + ph(0] p, /3)] _1 — d(pZ)/dp = 
(3aC e (a). 

The work presented here and elsewhere^ might have 
a didactic value in showing that relatively simple math- 
ematics allows one to check the thermodynamic self- 
consistency between the virial and energy routes in some 
important special cases. 
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